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RADIO IDENTIFICATION OF THE X-RAY JET IN THE Z=4.3 QUASAR GB 1508+5714 
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ABSTRACT 

The recent discovery of an X-ray jet in the z=4.3 quasar GB 1508+5714 by Yuan et al. (astro-ph/0309318) 
and Siemiginowska et al. (astro-ph/03 10241) prompted a search for its radio counterpart. Here, we report the 
successful discovery of faint radio emission from the jet at 1 .4 GHz using archival VLA data. The X-ray emission 
is best interpreted as inverse Compton (IC) emission off the CMB as discussed by the previous investigators. In 
this scenario, its high X-ray to radio monochromatic luminosity ratio, compared to previously detected IC/CMB 
X-ray jets at lower redshift, is a natural consequence of its high redshift. 

Subject headings: Galaxies: active — galaxies: jets — quasars: general — quasars: individual (GB 1508+5714) 
— radio continuum: galaxies — X-rays: galaxies 



1. BACKGROUND 

Since its launch in 1999, the Chandra X-ray Observatory 
has been used to detect a large number of X-ray jets in Active 
Galactic Nuclei (AGN), where prominent radio jets were pre- 
viously known to exist (see e.g., Harris & Krawczynski 2002, 
and associated website 1 ). The recent report (Yuan et al. 2003; 
Siemiginowska et al. 2003b) of an extended X-ray jet originat- 
ing from the z=4.3 quasar GB 1508+5714, where previous ob- 
servations showed no obvious sign of extended radio emission, 
presents an interesting case. The X-ray feature is strong - well 
over 100 counts were detected from it in the ^90 ksec Chandra 
exposure. An archival HST image helps rule out the possibility 
that it is due to a foreground galaxy or a gravitationally lensed 
image of the quasar (Siemiginowska et al. 2003b). Based on 
deep X-ray source counts, it has a low probability of being a 
random unassociated X-ray field source. 

As discussed by the previous authors, such detections of X- 
ray jets at large redshifts are actually to be expected as a natural 
consequence of the inverse Compton (IC) off the CMB model 
(e.g., Tavecchio et al. 2000; Celotti, Ghisellini, & Chiaberge 
2001). This is because the (1+z) 4 dependence of the CMB en- 
ergy density compensates for cosmological dimming of radia- 
tion, so that IC/CMB X-ray jets should remain detectable out 
to large cosmological distances (Schwartz 2002a). The model 
has been successfully applied to account for X-ray jets in many 
other powerful quasars at more modest redshifts (e.g., Sam- 
bruna et al. 2002), requiring that the jets are still highly rela- 
tivistic on kilo-parsec scales, in order that the electrons in the 
jet frame see an adequately boosted photon source. However, 
the lack of a detection of the GB 1508+5714 jet at lower fre- 
quencies, along with only a rough constraint on the X-ray spec- 
trum, could not rule out a synchrotron origin for the X-rays 
(Siemiginowska et al. 2003b). A previous search in the radio for 
a proposed X-ray jet in another high redshift (z=5.99) quasar, 
SDSS 1306+0356 (Schwartz 2002b), yielded only an upper 
limit of <0.1 mJy at 1.4 GHz (Petric et al. 2003; Schwartz, 
Cheung, & Wardle 2003). 

Distinguishing between the two possible emission processes 
is important, as they probe different energetic phenomena. In 



the case of synchrotron X-ray emission, the X-rays mark sites 
of particle acceleration with electrons accelerated up to 7 ~ 10 7 
with very short lifetimes (e.g. M87; Harris et al. 2003). When 
considered together with minimum energy /equipartition condi- 
tions, the IC/CMB model offers important constraints on the 
beaming, magnetic field, and jet power (Tavecchio et al. 2000; 
Celotti, Ghisellini, & Chiaberge 2001). 

Both the IC and synchrotron interpretations of X-ray jet 
emission usually require a population of relativistic electrons 
emitting synchrotron radiation at lower (radio) frequencies. 
Based on the X-ray flux and spectrum, the two models give 
different predictions of the radio component flux and spectrum. 
Here, we report the detection of such a radio component coin- 
cident with the X-ray feature extending from GB 1508+5714 
from an analysis of archival VLA data. This supports its in- 
terpretation as an X-ray jet and we discuss the X-ray emis- 
sion in the context of the new 1.4 GHz detection, along with 
the previously set optical, and a new 8.4 GHz limit. Fol- 
lowing Yuan et al. (2003) and Siemiginowska et al. (2003b), 
H = 71 km s" 1 Mpc" 1 , fJ M = 0.27 and n vac = 0.73 (Spergel et 
al. 2003) are assumed throughout, so 1" = 6.871 kpc. 

2. ARCHIVAL RADIO OBSERVATIONS 

Two datasets which utilized the VLA in its highest resolu- 
tion A-configuration were obtained from the NRAO 2 archive. 
The data were calibrated in the NRAO AIPS package (Bridle & 
Greisen 1994) and brought into DIFMAP (Shepherd, Pearson, 
& Taylor 1994) for editing, self-calibration, and imaging. The 
flux density scale was set on the VLA 1999.2 scale using scans 
of 3C 286 as outlined in the VLA Calibrator Manual (Perley & 
Taylor 2003). 

The 1 .4 GHz dataset - a single five minute snapshot obser- 
vation - yielded a 10cr detection of a 1.2 mJy feature extended 
from the quasar (Figure 1). This feature is only about 0.5% of 
the core flux (224 ±5% mJy), so it is no surprise that the orig- 
inal investigators (Moran & Helfand 1997) simply concluded 
that the quasar was an unresolved point source from these data. 
The measured off-source RMS of ,$,0.12 mJy in the image is 
within 50% of the thermal noise limit expected from the in- 
tegration time. The resultant image has a dynamic range of 



1 http://hea-www.harvard.edu/XJET/ 

2 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
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almost 2,000:1. 

The jet knot was modeled with an elliptical gaussian profile 
in both the (u,v) plane with DIFMAP's MODELFIT utility, and 
in the map plane with the JMFIT task in AIPS. The fits gave 
consistent measures of the deconvolved size: l"x < 1", elon- 
gated along the jet direction. We note that this is comparable 
to the beamsize for this observation so the source is consistent 
with being unresolved, and the size should be considered an 
upper limit. The separation from the core is 2.5" at a position 
angle of -114°, and is comparable to that measured from the 
X-ray data. Interestingly however, the peak in the X-ray jet ap- 
pears somewhat closer to the nucleus than the radio centroid in 
the image overlay (Figure 2), by about one Chandra ACIS pixel 
(0.492" per pixel), or equivalently ^3 kpc projected distance. 
The images were aligned by the cores to better than half of a 
pixel. Siemiginowska et al. (2003b) found the peak in the X-ray 
jet to be ^2" away (see also their subpixel re-binned radial pro- 
file published in their Figure 2), consistent with our qualitative 
assessment. Yuan et al. (2003) stated the same peak to be ^3" 
distant from the nucleus, although we measured a smaller value 
off their Figure 1 . Much more apparent X-ray /radio offsets are 
seen in other powerful quasars (e.g. up to ~2" in PKS 1127- 
145; Siemiginowska et al. 2002). Most, if not all of the X-ray 
counts from the GB 1508+5714 jet lie within the outermost ra- 
dio contour in the VLA image (Figure 2). 

A 10-minute 8.4 GHz observation obtained on Nov 3, 1996 
(Program AD388) yielded no detection of the radio jet. The 
measured off-source RMS in the naturally weighted image 
(^0.35" beam) is at about the expected thermal noise limit of 
0.045 mJy. In order to judge if the non-detection was a result 
of the greater resolution in this image compared to the 1 .4 GHz 
map, the 8.4 GHz dataset were tapered by different amounts. 
No outstanding feature appeared above the residual artifacts 
from the dirty beam. The measured RMS at the expected posi- 
tion of the jet in a tapered image restored with a 0.75" beam is 
0.3 mJy (3a). This, along with the 1.4 GHz detection, is consis- 
tent with an a ^ 0.8 radio spectrum (F v oc v~ a ), which agrees 
with the measured X-ray spectral index: 0.9±0.36 (Siemigi- 
nowska et al. 2003b) and 0.92^g;|| (Yuan et al. 2003). 

3. CONSTRAINTS FROM THE RADIO DATA 

The spectral energy distribution of the jet knot in 
GB 1508+5714 is shown in Figure 3. If the X-ray and ra- 
dio emission are drawn from the same population of relativis- 
tic electrons emitting synchrotron radiation, the spectral index 
will be about what is measured between the 1 .4 GHz and X-ray 
(using 1 .68 x 10~ 6 ph/cm 2 /s/keV at 1 keV, reported by Siemigi- 
nowska et al. (2003b)) detections: a rx = 0.73. This is a typi- 
cal value for the spectrum seen in radio jets (Bridle & Perley 
1984), and is consistent with the measured X-ray spectrum of 
~ 0.9 ±0.36 (Yuan et al. 2003; Siemiginowska et al. 2003b), 
and our constraint on the radio spectrum of ^0.8 (the mod- 
est 8.4 GHz dataset does not preclude that the higher frequency 
radio emission, >40 GHz in the source frame, was actually re- 
solved out). The optical HST limit (3cr) from Siemiginowska et 
al. (2003b) is not useful in this case, as it hovers over the radio- 
to-X-ray spectrum, and only weaker constraints are available at 
other optical bands (Yuan et al. 2003). 

We can estimate an equipartition magnetic field of B eq ~ 
9.7 x 10" 5 5 _1 ~f-° in n (r/ 1 ff 21 G, by adopting the observed a rx 
value as the optically thin spectral index, and that the particles 
and field fill a fraction / of a 2.9 x 10 67 cm 3 sphere, correspond- 



ing to the size upper limit derived from the radio detection. The 
factor r\ is the ratio of magnetic field and electron energy densi- 
ties, and 7 ml „ is the lower energy cutoff of a power law distribu- 
tion of relativistic particles. If we further assume r/=l, f=l, and 
7mm=10, then B eq ~ 7.4 x 10~ 5 <5 _1 G, where 8 is the unknown 
jet Doppler beaming factor. In this field, electrons emitting the 
highest energy X-rays detected (~5 keV; Siemiginowska et al. 
2003b) will have 7 <~ 10 8 and short lifetimes (^50 years) - this 
would require continuous in-situ reacceleration of high energy 
particles in order that the X-ray jet not be a transient feature. At 
this high redshift, the energy density of the CMB already ex- 
ceeds that of the equipartition field, even without bulk motion 
(see below), so inverse Compton losses will already be domi- 
nant. The lifetime of the 7 ~ 10 8 electrons calculated above is 
then a strict upper limit. Taking this evidence together, it is un- 
likely that the jet X-ray emission is dominated by synchrotron 
losses. 

An IC/CMB origin for this high redshift X-ray jet is preferred 
(Yuan et al. 2003; Siemiginowska et al. 2003b), and give us an 
additional constraint on the allowed range of magnetic field and 
jet Doppler factor to the equipartition condition. The following 
expression is taken from Tavecchio (2002): 



B = S[ 



3 c(o) 1 v" fi("c) -,-i/(i-ta) 



2 a T c U rad v%- 1 v? F s (v s ) 



(1) 



where, c(a) is a dimensionless function of a (e.g., Ghisellini, 
Maraschi, & Treves 1985), <ij is the Thomson cross section, 
and Urad = 4.19 x 10" 13 (1 +z) 4 ergs cm" 3 for the CMB (e.g., 
Schwartz 2002a), whose spectrum is peaked near vq = 1 .6 x 
10 11 (1+z) Hz (Harris & Krawczynski 2002). Assuming a = 1, 
the observed ratio of the Compton (X-ray) and synchrotron (ra- 
dio) fluxes (F c and F s , observed at v c and v s , respectively), allow 
us to estimate B ~ 7.2 x 10" 6 8 G. The corresponding equiparti- 
tion calculation gives B eq ~ 10" 4 8' 1 G, making the assumptions 
as above about i], f, and -f min . Both estimates of B-field are rea- 
sonably insensitive to the assumed spectral index: up to ~25% 
smaller for a = 0.9, the measured X-ray spectral index. In order 
to reconcile the differences between the two equations (oppo- 
site dependencies of B on 8), a Doppler factor greater than one 
is required. These two constraints bracket a plausible range of 
solutions around a 30^iG field and S ~4. Siemiginowska et al. 
(2003b) arrived at similar conditions using slightly different as- 
sumptions. The unresolved nucleus is a variable X-ray and (flat 
spectrum) radio source with a very high X-ray to optical lumi- 
nosity (see, Moran & Helfand 1997, and references therein). 
These observations have been taken as evidence of beaming on 
smaller scales, in the nucleus. 

The ratio of the monochromatic X-ray to radio luminosities 
(i.e., in vF v ) of the GB 1508+5714 jet is 158. This is one of the 
highest amongst the known X-ray jets which can be attributed 
to IC/CMB emission so far. In Figure 4, this ratio is plotted vs. 
redshift along with IC/CMB X-ray jets taken from the literature. 
The values seem to vary widely over 4 orders of magnitude, 
but predominantly show X-ray to radio ratios greater than 1. 
The jet in GB 1508+5714 sees 1-2 orders of magnitude times 
greater energy density from the CMB than other jets at lower 
redshift, so it is tempting to speculate that its extreme redshift 
may account, to first order, for its large X-ray to radio lumi- 
nosity ratio. This could similarly, account for the large X-ray 
to radio ratio limit obtained for the proposed X-ray jet in the 
z=5.99 quasar SDSS 1306+0356 (Schwartz 2002b; Schwartz, 
Cheung, & Wardle 2003, Fig. 4). 
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The observed monochromatic X-ray to radio luminosity ra- 
tios can be compared to what is expected in the IC/CMB model. 
We can obtain an expression for the expected ratio by rearrang- 
ing equation 1 (and setting a = 1): 

^ C J F> c ) = 4.19xlO- 13 (l+z) 4 
v, F s (v s ) (B/S) 2 /$ir ' 

The flux ratio is simply proportional to the ratio of the energy 
densities in the CMB and magnetic field, with a Doppler factor 
term. This expression is plotted in Figure 4 for several different 
combinations of B and S, one of which, best suits the observa- 
tions of GB 1508+5714. It appears that varying B and S can 
account for the large spread of observed X-ray to radio lumi- 
nosity ratios within a given redshift range. Larger B-field, or 
more likely, smaller S may account for the low values observed 
in Q0957+561 (z=1.41, Chartas et al. 2002) and 3C 9 (z=2.012, 
Fabian, Celotti, & Johnstone 2003), as they are known to have 
low radio core dominance, and weak VLBI structures (Camp- 
bell et al. 1995; Hough et al. 2002) compared to the other 
sources studied here. In all but 3C 179 (z=0.846, Sambruna et 
al. 2002) where multiple knot regions along the jet can be dis- 
tinguished, the X-ray to radio flux ratio decreases with increas- 



ing distance from the nucleus, and can similarly be accounted 
for by varying B/S. As discussed specifically in the case of 
3C 273 (z=0.158, Sambruna et al. 2001), the variations along 
the jet may indicate deceleration, or increasing B along the jet, 
possibly by compression in strong shocks. Future detections in 
the redshift range ~2-4 or higher with large X-ray/radio flux ra- 
tios can lend further support to the currently preferred IC/CMB 
model used to explain the X-ray emission in quasar jets. 
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FIG. 1 . — Naturally weighted VLA 1.4 GHz archival snapshot image of GB 1508+5714 (at center) and its ~2.5" long radio jet extending south of west of the 
quasar. The contours begin at 0.25 mJy/beam (2 times the measured RMS in the image) and the positive values (solid contours) are spaced by factors of V2, with 
an image peak of 224 mJy/beam. The restoring beam is plotted at the bottom left corner and has dimensions 1.8" X 1.35" at a position angle of -3.61°. 
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FIG. 2. — Chandra X-ray image of GB 1508+5714 (color) from the archive data published in Yuan et al. (2003) and Siemiginowska et al. (2003b), with VLA 1.4 
GHz image overlaid. The quasar is at the origin, and the jet feature indicated. The radio image is from the same data presented in Figure 1 but restored with the 
uniformly weighted beam (1.52" X 1.03" at PA=^t.44°). The lowest contour plotted is 0.2 mJy/beam and subsequently spaced by factors of two. 
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FIG. 3. — Radio-to-X-ray spectral energy distribution of the GB 1508+5714 jet. The bowtie indicates the X-ray spectrum measurement with Icr uncertainty 
from Siemiginowska et al. (2003b) and the arrows indicate 3<r limits. The radio data are from this work, and the X-ray point and optical limit were taken from 
Siemiginowska et al. (2003b). A dotted line shows the a rx slope and simply joins the radio and X-ray detections. 



Radio Jet in z=4.3 Quasar 




Redshift 



FIG. 4. — Plot of the ratio of the jet X-ray to radio monochromatic luminosity vs. redshift. Only jet features interpreted by the authors as IC/CMB X-ray emission 
are plotted. The curves indicate the expected ratio for given combinations of B and 8, which scale as (1+z) 4 . For reference, the B ~ 30/iG and <5 ~4 case derived for 
GB 1508+5714, which used the additional equipartition constraint, defines the dotted line which lies in between the other two curves. Light vertical lines connect 
features from the same source. The data are for (in order of increasing redshift) 3C 273-A to D (Sambruna et al. 2001), 1 150+497-A, B, PKS 1 136-135-B (Sambruna 
et al. 2002), B2 0738+3 13-A (Siemiginowska et al. 2003a), PKS 0637-752-WK7.8 (Chartas et al. 2000), 3C 207-knot (Brunetti et al. 2002), PKS 1354+195-A, 
B, 3C 179-A, B (Sambruna et al. 2002), PKS 1127-145-A to C (Siemiginowska et al. 2002), Q0957+561-B, C (Chartas et al. 2002), 3C 9-jet (Fabian, Celotti, & 
Johnstone 2003) with the corresponding radio flux (sum of knots F to J) taken from Bridle et al. (1994), and GB 1508+5714 and SDSS 1306+0356 (see text). 



